LIBRARY 

TECHNICAL    REPORT   SECTION 
NAVAL    POSTGRADUATE   SCHOOL' 
MONTEREY.  CALIFORNIA    93946 


NPS-67Gn7706J_ 


NAVAL  POSTGRADUATE  SCHOOL 

/y  Monterey,  California 


ON   THE   STABILITY  OF  TWO   BASIC   PARALLEL   FLOWS 

ANALYSIS,    PROGRESS    REPORT  AND  PROPOSAL 

FOR  FURTHER  RESEARCH 


T.    H.    GAWAIN 


JUNE   1977 


FEDDOCS 

D  208.14/2:NPS-67Gn77061 


NAVAL  POSTGRADUATE  SCHOOL 
Montere,  California 


Rear  Admiral  I.  W.  Linder  Jack  R.  Borsting 

Superintendent  Provost 


Reproduction  of  all  or  part  of  this  report  is  authorized. 


This  report  was  prepared  by: 


UNCLASSIFIED 


ECURlTY   CLASSIFICATION  OF   THIS  PACE  C*h—>  Omtm  gnl.f  d) 


REPORT  DOCUMENTATION  PAGE 


REPORT  NUMItft 

NPS-67Gn77061 


2.  GOVT   ACCESSION  NO. 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.     RECIPIENT'S  CATALOG  NUMBER 


.     TITLE  ^andSuMllft; 

ON  THE  STABILITY  OF  TWO  BASIC  PARALLEL  FLOWS 
ANALYSIS,  PROGRESS  REPORT  AND  PROPOSAL  FOR 
FURTHER  RESEARCH 


K     AUTHORS 

T.  H.  GAWAIN 


5.  TYPE  OF  REPORT  *  PERIOO  COVERED 


t.  PERFORMING  ORG.  REPORT  NUMBER 


«.  CONTRACT  OR  GRANT  NUMBERf*) 


r  PERFORMING  ORGANIZATION  NAME  AND  AODREJS 

Naval  Postgraduate  School 
Monterey,  California  93940 


X 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  a   WORK  UNIT  NUMBERS 


I.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 


12.     REPORT   DATE 

1  June  1977 


11.     NUM3E*  IP  PAGES 

61 


U.     MONITORING  AGENCY   NAME  *    AOORE5SC</  dltt.r~,t  trvm  Controlling  0»(c> 


IS.     IZCURIT-'   CLASS.  (ol  trU*  f*po/t) 


TTMPT  A^qTFTF.n 


lla      DECL  ASSIFICAT'.ON/ V,u*NGRAOlNG 
SCHEDULE 


16.     DISTRIBUTION   STATEMENT  (el  t*lt  Rrpott) 

Approved  for  public  release;  distribution  unlimited. 


,7,     DISTWfU       ON  STATEMENT  (ol  fh.  O.frtttf  -f~*  '»  *'•«*  20.  II  dHt—»»  *~>  *•*>«) 


18.     SUPPLEMENTARY  NOTES 


1».     KEY   WORDS  (Ctittmt*  an 


'~4~  II  n»cfTT  n*  Id^tlty  by  k/o«*  numkt) 


Hydrodynamic  Stability 
Plane  Poiseuille  Flow 
Pipe  Poiseuille  Flow 


20       ABSTRACT  (Centime  on  r»*~f  •  !*•  It  n.C««— T  -**  <*«•<'*  *7  «•«* 

This  report  provides  a  detailed  technical  outline  and  evaluation  of 
recent  research  performed  at  the  Naval  Postgraduate  School  on  the  hydrodynamic 
stability  of  plane  Poiseuille  flow  and  of  pipe  Poiseuille  flow.   These  two 
cases  are  of  fundamental  importance  in  connection  with  the  theory  of  hydro- 
dynamic  stability.   Each  case  involves  significant  discrepancies  between  the 
predictions  of  conventional  theory  and  the  results  of  actual  experimental 
observations.   In  particular,  the  conventional  theory  fails  completely  to 


FORM 


DD     IJJmTtI     1473  EDITION   OF    I   NOV  •»  IS  OBSOLETE 

(Page      1)  S/N   0102-014- 6601   I 


TTNrT.ASSTFTFTi 


SECURITY  CLASSIFICATION  OF   TH.S  PAGE  (**•*  D.i.  «»•« 


UNCLASSIFIED 


fliCUWTV   CLASSIFICATION   OF    THIS  P»GEfW>..n   O-tm  Enfr.i 


account  for  the  well  known  instability  of  ordinary  pipe  flow. 

This  report  describes  a  more  general  theory  which  shows  promise  of 
overcoming  the  above  limitations.  The  new  theory  involves  a  number  of 
innovations  in  the  formulation  of  some  of  the  basic  equations,  in  the  fon 
lation  of  certain  boundary  conditions,  and  in  the  formulation  of  the 
criterion  of  stability  itself.  The  analytical  work  is  essentially  complet 
and  preliminary  calculations,  while  still  of  limited  scope,  appear  to 
support  the  theory  and  are  therefore  quite  encouraging. 

Nevertheless,  before  the  new  theory  can  be  considered  as  definitive!] 
confirmed,  a  considerable  amount  of  systematic  computer  calculation  will 
still  be  necessary.   Research  support  is  now  being  sought  for  this  purpos 


DD     Form       1473 

.  1  Jan  73 
S/N    0102-014-6601 


UNCLASSIFIED 


SECURITY   CLASSIFICATION  OF  THIS  PAGEr***"  Dmt  I 


111 


SUMMARY 

This  report  provides  a  detailed  technical  outline  and  evaluation  of 
recent  research  performed  at  the  Naval  Postgraduate  School  on  the  hydro- 
dynamic  stability  of  plane  Poiseuille  flow  and  of  pipe  Poiseuille  flow. 
These  two  cases  are  of  fundamental  importance  in  connection  with  the 
theory  of  hydrodynamic  stability.  Each  case  involves  significant  discrep- 
ancies between  the  predictions  of  conventional  theory  and  the  results  of 
actual  experimental  observations.   In  particular,  the  conventional  theory 
fails  completely  to  account  for  the  well  known  instability  of  ordinary 
pipe  flow. 

This  report  describes  a  more  general  theory  which  shows  promise  of 
overcoming  the  above  limitations.   The  new  theory  involves  a  number  of 
innovations  in  the  formulation  of  some  of  the  basic  equations,  in  the  formu- 
lation of  certain  boundary  conditions,  and  in  the  formulation  of  the 
criterion  of  stability  itself.  The  analytical  work  is  essentially  complete 
and  preliminary  calculations,  while  still  of  limited  scope,  appear  to 
support  the  theory  and  are  therefore  quite  encouraging. 

Nevertheless,  before  the  new  theory  can  be  considered  as  definitively 
confirmed,  a  considerable  amount  of  systematic  computer  calculation  will 
still  be  necessary.   Research  support  is  now  being  sought  for  this  purpose. 
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1.   Introduction 

The  basic  research  discussed  in  this  report  deals  with  the  classical 
problem  of  the  hydrodynamic  stability  of  parallel  flows.   Two  fundamental 
cases,  both  of  considerable  historical  and  practical  importance,  are  of 
particular  interest,  namely,  the  plane  Poiseuille  flow  between  infinite 
parallel  plates  and  the  flow  through  a  uniform  pipe.   Despite  the  consid- 
erable amount  of  research  that  has  been  done  over  the  years  on  these  two 
fundamental  cases,  the  theory  available  up  to  now  does  not  account  ade- 
quately for  the  observed  experimental  facts.   Indeed,  for  the  case  of  pipe 
flow,  the  classical  theory  fails  to  predict  the  onset  of  hydrodynamic 
instability  at  any  Reynolds  number  whatever!   In  this  respect  the  existing 
theory  is  clearly  inadequate. 

A  critical  review  of  the  situation  by  the  present  investigator  has  dis- 
closed that  the  above  limitations  and  defects  of  the  classical  theory  can 
largely  be  overcome  by  generalizing  the  analysis  to  include  certain  three 
dimensional  effects  that  are  neglected  in  the  classical  approach.   The  re- 
sulting calculations  have  succeeded  for  the  first  time  in  revealing  instabi- 
lities in  pipe  flow  for  certain  types  of  three  dimensional  perturbations; 
these  preliminary  results  appear  to  be  consistent,  at  least  qualitatively, 
with  the  known  experimental  facts  about  pipe  flow.   The  corresponding 
calculations  for  plane  Poiseuille  flow  reveal  that  under  certain  conditions 
instability  can  occur  at  a  Reynolds  number  significantly  lower  than  that 
predicted  by  the  classical  theory.   The  new  analysis  therefore  promises  to 
advance  our  basic  understanding  of  the  important  phenomena  of  hydrodynamic 
instability  in  a  significant  way. 
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The  exploration  of  the  quantitative  implications  of  this  generalized 
three  dimensional  theory  entails  a  heavy  burden  of  computer  calculations. 
Research  support  is  needed  for  the  purpose  of  carrying  out  systematic 
calculations  of  this  kind.   Support  is  also  needed  for  the  purpose  of 
documenting  and  publishing  in  the  technical  literature  the  analytical 
basis  of  the  theory  and  such  preliminary  numerical  results  as  have  already 
been  obtained. 

More  specifically,  it  is  proposed  to  carry  out  the  above  documentation 
and  publication,  as  well  as  a  significant  amount  of  the  needed  additional 
calculations  during  the  first  year  of  the  proposed  research.   If  warranted 
by  the  results  accomplished  during  the  first  year,  additional  support  for 
follow  on  research  will  probably  be  sought  also  for  another  year  or  two 
thereafter,  to  consolidate  the  progress  made. 
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2.   Analytical  Approach 

The  essentials  of  the  analytical  approach  involved  can  best  be  ex- 
plained by  first  considering  the  specific  case  of  plane  Poiseuille  flow. 

All  equations  are  nondimensionalized  by  utilizing  as  dimensional 
reference  parameters  the  density  p  of  the  fluid,  the  semi-height  a 
between  the  walls  and  the  mean  volumetric  velocity  U  of  the  fluid. 

Cartesian  coordinates  x,  y,  z  are  used  with  x  in  the  direction  of 
the  mean  flow  and  y  normal  to  the  walls.   Dimensionless  coordinate  y 
varies  from  -1  at  the  lower  wall  to  +1  at  the  upper  wall.   Symbols 
i,  j,  k  denote  unit  vectors  in  the  x,  y,  z  directions,  respectively. 

It  is  well  known  that  the  dimensionless  velocity  U  of  the  mean  flow 
along  the  x  axis  is  given  by  the  simple  expression 

U  =  |  (1  -  y2)  (2.1) 

The  hydrodynamic  stability  of  this  flow  field  can  be  analyzed  in  the 
usual  way  by  superimposing  upon  the  above  mean  flow  a  suitable  perturbation 
velocity  v  of  small  amplitude.   It  is  of  course  necessary  that  v  satisfy 
the  continuity  equation  for  incompressible  flow,  namely, 

V'v  =  0  (2.2) 

A  convenient  way  to  ensure  that  Eq.  (2.2)  is  satisfied  is  to  express 

v  in  terms  of  a  corresponding  vector  potential  function,  let  us  call  it 

■*■ 

W  ,  by  a  relation  of  the  form 

v  =  VXW  (2.3) 

It  is  customary  in  linear  stability  analysis  to  represent  quantities 
like  v  and  W  by  complex  variables  with  the  understanding  that  only 
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the  real  part  of  the  variable  shall  be  identified  with  the  actual  physical 
effect  in  question.   This  convention  simplifies  the  analytical  form  of 
the  equations  involved. 

For  a  given  perturbation  vector  field  v  ,  the  corresponding  vector 

-> 
potential  function  W  is  not  uniquely  determined  but  may  be  expressed  in 

a  number  of  alternative  forms  all  of  which  are  ultimately  equivalent.   In 

the  present  case,  however,  it  proves  advantageous  to  represent  the  complex 

vector  potential  function  W  as  follows 

W  =  [i(o)  +  J  G(y)  +  k  H(y)]  e  (aX  +  Bz  +  yt)  (2.4) 

where  the  complex  constants  a,  3,  Y  mSLY   De  further  expanded  to  the  form 

a  =  a  +  i  a 
K       J. 

3  =  3R  +  i  3X  (2.5) 

y  =  yr  +  i  Yx 

and  where   i  denotes  the  imaginary  unit,  that  is, 

i  =  +  v^l  (2.6) 

The  quantities   G(y)   and  H(y)   in  Eq .  (2.4)  amount  to  complex  stream 
functions.   Although  a  single  stream  function  suffices  to  represent  any 
incompressible  plane  flow,  two  stream  functions  are  required  to  represent 
an  arbitrary  incompressible  three  dimensional  flow  with  complete  generality. 

Incidentally,  this  use  of  two  stream  functions  to  represent  the  pertur- 
bation flow  in  this  problem  is  an  innovation  introduced  by  the  present 
investigator.   As  far  as  he  knows,  this  method  does  not  appear  in  the 
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existing  technical  literature  which  deals  with  the  present  problem. 

Of  course,  the  actual  forms  of  the  complex  functions  G(y)   and  H(y) 
are  initially  unknown  but  are  ultimately  calculable  in  the  course  of  the 
final  numerical  solution. 

The  notation  used  for  the  exponential  coefficients  in  Eqs.  (2.4)  and 
(2.5)  was  adopted  for  reasons  of  overall  generality,  consistency  and 
simplicity  in  the  subsequent  analysis.   Various  alternative  but  equivalent 
notations  are  also  found  in  the  technical  literature,  but  all  such  alter- 
natives can  be  readily  translated  to  the  above  format  whenever  necessary. 

Parameters  a   and  S   in  Eq.  (2.5)  are  associated  with  the 
imaginary  unit  i  ,  and  hence  they  represent  the  oscillatory  wave  numbers 
of  the  perturbations  with  respect  to  space  coordinates  x  and  z,  re- 
spectively.  On  the  other  hand   a   and   6^  do  not  involve  i  and  hence 

R         R 

they  denote  real  exponential  growth  rates  (or,  if  negative,  decay  rates) 

with  respect  to  space  coordinates  x  and  z  .  Likewise,  y       represents 

the  oscillatory  wave  number  of  the  perturbation  with  respect  to  time  t  , 

that  is,  the  angular  frequency,  while  Yt>  represents  the  exponential 

growth  or  decay  rate  with  respect  to  time.   Incidentally,  it  is  convenient 

to  restrict  a   and   BT   to  non-negative  values  and  there  is  no  loss  of 

generality  involved  in  doing  this  so  long  as  a   and  Q       remain  free  to 

R         R 

take  on  positive  or  negative  values  as  required. 

Notice  that  the  above  quantities  and  equations  are  defined  initially 

in  relation  to  axes  that  remain  fixed  with  respect  to  the  walls.   Later  we 

shall  also  consider  axes  that  translate  at  constant  velocity  relative  to 

the  fixed  walls.   Such  motion  does  not  affect  parameters  a  ,  a  ,  BD ,  8 

R    I    R    I 

but  can  in  general  change  the  apparent  values  of  y        and  y      .   We  denote 

R         1 

these  changed  values  with  respect  to  the  moving  axes  by  symbols  y'   and  y' 

R  1 
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In  the  method  of  analysis  considered  here,  the  numerical  values  of 

a  ,  a  ,  8  ,  8T   are  assigned  arbitrarily,  but  the  corresponding  values  of 
K    1    K    1 

Y   and  yT  must  be  found  by  calculation.   This  point  will  be  clarified 
K         J. 

further  in  the  subsequent  discussion. 

Of  course  it  is  essential  that  the  resultant  perturbed  flow  field 
defined  by  the  foregoing  relations  satisfy  Newton's  second  law  of  motion. 
This  in  turn  requires  that  the  velocity  and  vorticity  components  satisfy 
the  vorticity  transport  equation,  which  is  the  relation  that  is  found  by 
applying  the  curl  operator  to  the  vector  equation  of  motion.   This  relation 
involves  terms  which  are  linear  in  the  perturbation  quantities  and  terms 
which  are  quadratic  in  the  perturbation  quantities.   For  small  perturba- 
tions, however,  the  quadratic  terms  may  be  neglected.   In  this  way  the 
customary  linearized  equation  is  obtained.   A  full  explanation  of  these 
details  lies  outside  the  scope  of  the  present  discussion,  but  is  given  by 
Harrison,  Ref.  (6).   Fortunately,  the  result  is  well  known.   One  convenient 
version  of  this  result  may  be  summarized  concisely  as  follows 

~  VX  (Vxw)  +  Vx  [3  x  v  -  U  x  5]  +  (J^J   =0  (2.7) 

where  the  symbols  are  defined  as  indicated  below. 

U  =  ir  (1-y  )  Velocity  of  mean  flow  field 

u   =   VxU  Vorticity  of  mean  flow  field 

Perturbation  velocity  (2.8) 

Perturbation  vorticity 

Reynolds  number 
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V    = 

VxW 

-> 

OJ     = 

7xv 

R     ■ 

e 

Notice  that  Eq,  (2,7),  the  vorticity  transport  relation,  introduces 

yet  another  important  parameter  into  the  analysis,  namely,  the  Reynolds 

numb  er  R 
e 

Of  course  the  solution  of  Eq,  (2,7)  must  be  such  as  to  satisfy  the 
boundary  conditions  of  the  problem  which  are  simply  that 

v  =  0  at  y  =  +  1  (2,9) 

By  a  lengthy  but  systematic  procedure  which  is  fully  documented  in 
Ref  (6),  Eq.  (2.7)  can  be  expanded  in  terms  of  W  and  hence  in  terms  of  the 
functions  G(y)   and  H(y)  .   Eq.  (2,7)  is  a  vector  relation  and  therefore 
implies  three  corresponding  scalar  equations.   It  can  be  shown,  however, 
that  only  two  of  these  are  independent.   Consequently,  the  final  problem 
reduces  to  two  differential  equations  in  the  two  unknown  functions  G(y) 
and  H(y)  . 

In  this  connection,  the  following  three  auxiliary  functions  occur. 
These  functions  are  immediately  known  as  soon  as  definite  values  are 
assigned  to  parameter  R   and  to  the  complex  cons tans t  a  and   8  .   It 
is  helpful  to  calculate  these  functions  first,  before  proceeding  with  the 
rest  of  the  solution.   Thus 

T(y)  .  llL+jb.    .  „  1  (1-y2)   ■ 
e 

2    2 
Ity)  -  (a  |  g  }   +T(y)  (2.10) 

e 

J(y)  =  (a2  +  32).T(y)  -  3a 

Using  symbol   D   to  denote  the  differential  operator  —  ,  we  find  that 
the  resulting  vorticity  transport  equation  and  its  associated  boundary 
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condition  can  now  be  summarized  concisely  in  the  form 


j4-  D4H  +  I(y)  D2H  +  J(y)  h[  -  y   Jd2H  +  (a2  +  62)  H{  -  0 

where    H  =  0  \ 

I   at  y  =  +  1  (2.11) 

and     DH  =  0 


The  second  vorticity  transport  equation  and  its  associated  boundary 
condition  can  likewise  be  written 


j|-D2G  +  [T(y)  -  y]   g| 

D3H  +  [T(y)  -  y]  DH  -  3ayH[  (2.12) 


(a2  +  B2)  /  Re 


where   G  =  0   at  y  =  +  1 

Now  consider  the  situation  in  which  the  values  of  parameters 

a  ,  a  ,  3  ,  3T  and  R   are  arbitrarily  assigned.   Then  the  three  auxilary 
R    X    R    1         e 

functions  defined  by  Eqs .  (2.10)  become  known  functions.   Hence  the  basic 
relations  given  in  Eqs.  (2.11)  can  now  be  solved  numerically,  subject  to 
the  boundary  conditions  shown. 

Eqs.  (2.11)  define  an  eigenvalue  problem.   Theoretically  there  exist 
an  infinite  number  of  distinct  solutions  to  this  equation,  each  characterize 
by  a  corresponding  value  of  the  complex  constant  y    .   We  term  these  roots 
eigenvalues  and  denote  them  by  symbols  Ti»  To>  T*j  •••  •   Associated  with 
each  root  is  a  corresponding  complex  eigenf unction  H  (y) ,  H  (y) ,  H  (y) . . , 

Fortunately,  only  the  solution  which  corresponds  to  the  least  stable  root 
is  required  for  the  stability  analysis.  The  least  stable  root  is  defined 
in  more  detail  later. 
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Once  the  kth  root   y   and  the  kth  eigenfunction  H^(y)   are 
known,  Eqs.  (2.12)  may  be  solved  for  the  corresponding  complementary 
eigenfunction  G,  (y)  ,  subject  to  the  boundary  condition  shown. 

Notice  that  Eqs.  (2.11)  and  (2.12)  for  H   and  G  ,  respectively, 
are  uncoupled  and  may  therefore  be  solved  sequentially.   This  fact  simpli- 
fies the  solution.   Usually  we  are  interested  only  in  finding  the  least 
stable  root.   In  that  case  it  is  unnecessary  to  solve  Eqs.  (2.12). 

By  utilizing  the  method  of  finite  differences,  it  is  possible  to 
reduce  the  foregoing  equations  to  matrix  format.   The  resulting  matrix 
represents  a  standard  eigenvalue  problem  which  can  be  solved  numerically 
by  familiar  and  well  established  techniques. 

The  details  of  this  reduction  to  matrix  format  are  presented  in 

Ref.  (6)  and  are  not  repeated  here. 

The  real  part  of  each  eigenvalue  is  denoted  by  y   and  the  imaginary 

R 

part  by  y   .   For  reasons  that  are  more  fully  explained  later,  the  eigen- 
value that  has  the  algebraically  highest  value  of  y   is  termed  the  least 

R 

stable  root.   We  denote  this  highest  real  value  by  symbol  y   .   Of  course 

RL 

y    depends  on  the  particular  numerical  values  assigned  to  parameters 
RL 

a_,  a_.  B„.  BT   and  R   .  We  may  symbolize  this  relation  in  the  form 
R   I   R   I        e 

YRL  =  YRL  !V  "l-  SR>  6I>  V  (2'13) 

The  stability  of  the  flow  field  is  determined  by  the  numerical  value 

of  y    according  to  the  criteria  considered  in  the  next  two  sections. 

RL 
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3.   Conventional  Criterion  of  Stability 

Historically,  the  analysis  of  stability  has  concentrated  attention 
chiefly  on  the  important  but  nevertheless  special  case  in  which  the  spatial 
form  of  the  pertubations  is  taken  to  be  purely  sinusoidal.  This  amounts 
to  restricting  the  analysis  to  the  particular  case  for  which,  in  the  present 
notation, 

aR2  +  3R2  =  0  (3.1) 

Moreover,  any  time  that  Eq.  (3.1)  applies,  the  well  known  theorem  of 
Squire,  Ref.  (1),  also  applies.   When  expressed  in  terms  of  the  present 
nomenclature,  Squire's  theorem  simply  asserts  that  for  specified  values  of 
parameters  R   and  a   ,  the  most  unstable  condition  will  always  be  the 
one  for  which 

6I  =  0  (3.2) 

Thus  of  the  five  parameters   R  a  ,  a  ,  6  ,  Q       which  occur  in  the 

e  R   I   R   l 

general  case,  the  constraints  expressed  by  Eqs.  (3.1)  and  (3.2)  eliminate 
three.   Therefore  the  conventional  analysis  finally  retains  only  two  para- 
meters, namely,   R   and  a_  . 

e        I 

In  the  conventional  nomenclature  the  quantity  we  have  called  aT   is 
usually  represented  by  the  unsubscripted  symbol  a  ;  to  avoid  ambiguity  in 
the  present  more  general  discussion,  however,  it  is  important  for  clarity  to 
adhere  consistently  to  the  subscripted  notation  as  introduced  originally  in 
connection  with  Eqs.  (2.5). 

Under  conditions  where  Eq.  (3.1)  applies,  that  is,  for  the  conventional 
analysis,  the  stability  of  the  flow  field  is  determined  by  a  very  simple 
rule.   The  flow  is  stable,  or  at  least  neutrally  stable,  provided  that 
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If  the  inequality  applies,  the  perturbation  amplitude  decays  expo— 
nentially  over  time  and  the  flow  is  stable.   If  the  equality  applies,  the 
perturbation  amplitude  neither  grows  nor  decays  overtime  and  the  flow  is 
neutrally  stable.   If  neither  of  the  above  conditions  is   satisfied,  the 
perturbation  amplitude  grows  with  time  and  the  flow  is  unstable. 

Numerical  calculations  based  on  the  conventional  restriction  of  Eq, 
(3.1),  on  the  application  of  Squire's  theorem  as  expressed  by  Eq.  (3.2), 
and  on  the  simple  stability  rule  summarized  in  Eq.  (3.3),  permit  the  flow 
stability  to  be  determined  for  arbitrarily  assigned  values  of  R   and 
a   .   The  region  of  instability  in  the  R   vs  a   plane  has  the  general 
form  sketched  qualitatively  in  Fig.  3.1.   Notice  that  there  exists  a 
critical  Reynolds  number  below  which  sinusoidal  streamwise  perturbations 
of  all  wave  numbers  a   are  stable.   The  conventional  analysis  gives  the 
value  of  the  critical  Reynolds  number,  based  on  channel  semi-height  and 
volumetric  mean  velocity  as  about  3850,  Thomas,  Ref .  (2).   Of  course,  all 
of  these  conventional  results  are  familiar;  they  have  been  known  for  some 
time. 
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\.      Generalized  Criterion  of  Stability 

The  simple  stability  rule  summarized  in  Eq.  (3.3)  of  the  previous 
section,  although  correct  for  the  special  case  of  purely  sinusoidal 
Dscillations,  is  not  adequate  for  dealing  with  the  more  general  situation 
tfhere 

(<xR2  +  6R2)  #0  (4.1) 

Moreover,  under  these  conditions,  Squire's  theorem  no  longer  holds 
either,  so  that 

gj  t   0  (4.2) 

Thus  of  the  five  parameters  R  ,  o.       aT,  g_,  Sx  which  occur  in  the 

e   K   l   R   I 

jeneral  case,  none  can  now  be  dropped.   Compared  with  the  conventional 
:ase  which  involves  only  the  two  parameters  R   and  a   ,  this  represents 
i  considerable  increase  in  complexity.  The  additional  parameters  which 
>ccur  in  the  generalized  problem  are  associated  with  certain  instabilities 
?hich  are  not  revealed  in  the  conventional  analysis.  This  is  the  basic 
reason  why  the  generalized  analysis  gives  results  in  better  agreement 
/ith  experiment  than  does  the  conventional  method. 

While  the  simple  criterion  of  stability  expressed  by  Eq.  (3.3)  is 
iniversally  accepted  as  applying  to  the  conventional  case,  there  is  no 
such  unanimity  among  different  investigators  concerning  the  appropriate 
stability  criterion  to  apply  to  the  generalized  problem. 

The  analysis  of  the  previous  section  has  shown  that  if  we  arbitrarily 
Impose  the  restrictions  that 

aR=  0 

6R  =  0  (4.3) 

YRL  =  0 
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these  conditions  are  sufficient  to  ensure  neutral  stability.   Eqs.  (4.3) 
are  frequently  misinterpreted  as  representing  the  conditions  which  are 
necessary  to  ensure  neutral  stability.   In  fact  Eqs.  (4.3),  while  un- 
questionably sufficient,  are  not  necessary  for  this  purpose.   If  they 
were  necessary  conditions,  they  would  of  course  rule  out  of  consideration 
altogether  the  generalized  case  which  is  the  very  subject  of  the  research 
under  discussion.  Nevertheless,  this  seems  to  be  a  fairly  common  mis- 
understanding of  the  situation,  a  fact  which  might  help  explain  why  the 
fully  general  three  dimensional  case  has  not  received  the  degree  of 
research  attention  that  it  deserves. 

A  useful  clue  to  a  more  satisfactory  formulation  of  a  generalized 
stability  criterion  is  provided  by  the  existance  of  two  special  cases  in 
which  there  seems  to  be  general  agreement  on  how  the  stability  should 
properly  be  defined,  Refs.  (4)  and  (5).   These  are  summarized  in  Eqs. 
(4.4)  and  (4.5)  below. 

The  flow  field  is  stable,  or  at  least  neutrally  stable, 


if       aR  =  0 

when    y   <  0 
RL  — 


or  if    y     =  0 
R 

when    a     0 


(4.4) 


(4.5) 


A  significant  feature  of  Eqs.  (4.4)  and  (4.5)  is  that  they  both  refer 

explicitly  to  parameter  a   but  not  to   3   .  Thus  stability  is  defined 

K  R 

in  terms  of  the  component  of  spatial  growth  rate  in  the  streamwise  direction 
but  not  in  terms  of  the  component  in  the  transverse  direction  (if  any) . 
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In  analytical  work  such  as  considered  primarily  in  this  discussion, 
parameter  a  can  be  set  to  zero  independently  as  required  by  the  first 
of  Eqs.  (4,4)  and  8       and  the  other  parameters  can  be  chosen  arbitrarily. 

RL  for  the  least  stable  root  can  be  found  by  calculation, 

and  the  resulting  stability  determined  from  the  second  of  Eqs,  (4.4). 

While  this  method  is  favored  by  some  investigators,  it  fails  for  the 

important  case  where  a  4   0  . 

R 

In  a  hypothetical  experiment  involving  steady  oscillations  we  would 

have  y  =  0  in  accordance  with  the  first  of  Eqs.  (4.5).   The  observed 
R 

streamwise  growth  or  decay  rate  would  then  fix  the  value  of  a    for  the 

least  stable  root.   The  second  of  Eqs.  (4.5)  would  then  fix  the  stability. 

While  this  concept  can  hardly  be  faulted  for  the  special  case  to  which  it 

applies,  unfortunately  it  fails  for  the  important  case  where  yR  4   0  . 

One  basic  reason  for  the  diversity  of  stability  criteria  utilized  by 

different  researchers  is  that  in  the  general  case  the  apparent  exponential 

growth  or  decay  rate  of  the  perturbation  with  respect  to  time,  as  denoted 

by  parameter  y'  ,  depends  in  part  on  the  motion,  if  any,  of  the  observer's 
K 

reference  frame.   To  see  this  consider  axes  that  translate  at  constant 

velocity  parallel  to  the  fixed  walls .   Let  x   denote  the  translation 

velocity  component  in  the  downstream  or  x  direction,  and  let  z 

denote  the  velocity  component,  if  any,  in  the  lateral  or  z  direction. 

Let  xT,  y,  z'   denote  the  coordinates  of  a  point  fixed  with  respect  to  the 

moving  axes,  and  let  x,  y,  z  denote  the  coordinates  of  the  same  point 

with  respect  to  fixed  axes,   Then 

x  =  x'  +  X   t 

(4.6) 

z  =  z'  +  z   t 

o 

4-3 


The  exponent  in  Eq.  (2.4)  can  now  be  written  as  follows 
X  =  ax  +  3z  +  yt 


=  a(xf  +  x  t)  +  6(z'  +  z  t)  +  yt 
o  o 


ax'  +  3z'  +  (y  +  ax  +  Bz  )t 

o     o 


=  axT  +  8z'  +  Y't 


(4.7) 


where 


Y'  =  Y  +  ax  +  6z  (4.8) 

o     o 


Breaking  this  into  real  and  imaginary  parts  gives 


YR  =   YR+aRXo   +   6RZo 


Yi  =  YI  +  ai  Xo  +  6I   Zo 


(4.9) 


Thus,  with  respect  to  a  reference  frame  which  moves  with  any  specified 

values  of  x   and  z   ,  the  flow  field  may  be  said  to  be  stable,  or  at 
o        o 

least  neutrally  stable,  provided  that  for  the  least  stable  root 

The  least  stable  root  is  that  root  which,  for  fixed  values  of  the 

other  parameters,  has  the  algebraically  largest  value  of  y1      as  defined 

R 

above. 

Incidentally,  for  purely  sinusoidal  oscillations,  that  is,  for 


. 


Eq.    (4.10)    reduces   to 


aR=   6R=  0  (4.1i: 


Yli  =   Y^     °  <4'12 


4-4 


In  other  words  the  stability  is  the  same  with  respect  to  all  uniformly 
translating  reference  frames  in  this  case,  irrespective  of  the  values  of 
x   and  z   involved.   Thus  the  criterion  of  stability  is  simple  and 
unambiguous  under  these  circumstances. 

When  we  revert  to  the  general  case,  however,  the  apparent  stability 


•    • 


as  expressed  by  Eq.  (4.10)  depends  on  the  velocities  x  ,  z   of  the 
chosen  reference  axes.   Thus  the  concept  of  instability  is  not  completely 
and  uniquely  defined  until  parameters  x   and  z   are  designated  in  an 
appropriate  way. 

In  the  present  context  it  is  reasonable  to  require  that  these  velo- 
cities be  specified  in  such  a  way  that  the  resulting  stability  criterion 
be  consistent  with  Eqs .  (4.4)  and  (4.5)  for  the  two  special  cases  to  which 
these  relations  apply.   Comparison  of  Eq.  (4.10)  with  Eqs.  (4.4)  and  (4.5) 
reveals  that  this  will  indeed  be  the  case  provided  that  we  set 


z   =  0 
o 


x  >  0 
o 


(4.13) 


Consequently,  Eqs.  (4.9)  and  (4.10)  now  simplify  to  the  form 

Y*i  =  YRL  +  aR  *o  -  ° 


^iL   =  YIL  +  al  Xo 


(4.14) 


where  the  magnitude  of  x   still  remains  to  be  defined.   In  this  connec- 

* 

tion,  a  method  that  is  followed  by  some  investigators  is  to  choose  x 

o 


in  such  a  way  as  to  give 


Vzl   =  0  (4.15) 
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The  corresponding  value  of  x   ,  which,  is  now  termed  the  phase  velocity. 


is  seen  to  be 


x  = 

o      aT 


(4.16) 


and  this  quantity  is  normally  positive  as  required.   Consequently,  with 
respect  to  this  particular  reference  frame,  the  flow  is  stable,  or  at  least 
neutrally  stable,  provided  that 

It  is  also  possible  to  define  a  phase  velocity  which  satisfies  Eq, 
(4.15)  in  a  more  general  manner.   In  place  of  setting   z  =  0  ,  we  may  re- 


quire that 


.  2    •  2 

x   +  z    =  a  minimum  (4.18) 

o     o 


However,  this  option  contravenes  Eq.  (4.13)  and  will  therefore  not  be  con- 
sidered further. 

Despite  the  fact  that  Eq.  (4.17)  is  the  stability  criterion  preferred 
by  some  writers,  this  investigator  does  not  recommend  its  use,  for  the 
following  reason.   Note  that  the  restriction  imposed  by  Eq.  (4.15)  is 
essentially  arbitrary.   It  is  a  restriction  on  the  angular  frequency  and 
as  such  has  no  necessary  relation  to  growth  or  decay  rate.   Consequently 
the  stability  criterion  expressed  by  Eq,  (4. 17), while  mathematically  correct 
with  respect  to  the  axes  to  which  it  pertains,  is  as  arbitrary  as  are 
those  axes  themselves. 

To  eliminate  the  arbitrary  element  involved  in  Eq.  (4.17),  this  investi- 
gator advocates  a  criterion  based  on  adopting  a  Lagrangian  viewpoint,   Con- 
sider a  given  streamline  of  the  mean  flow  located  at  some  specified  value 
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pf  coordinate  y  .   The  fluid  particles  which,  lie  along  this  streamline 
advance  downstream  at  an  average  velocity  U  relative  to  fixed  axes. 
We  now  stipulate  that  the  flow  along  this  streamline  shall  be  deemed 
stable,  neutrally  stable  or  unstable,  respectively,  according  to  whether 
the  particles  advancing  along  that  streamline  undergo  gradually  in- 
creasing,  unchanging  or  gradually  decreasing  amplitudes  of  oscillation. 
This  amounts  to  letting  the  references  axes  move  with  the  mean  velocity 
of  the  fluid  particles  along  the  chosen  streamline.  Hence  we  set 

xn  =  U(y)  =  |  (1-y2) 

(4.19) 

z   =  0 
o 

Consequently,  the  flow  along  the  given  streamline  maybe  said  to  be 
stable,  or  at  least  neutrally  stable,  provided  that  for  the  least  stable 
root 

TRL  =  YRL  +  aR  U(y)  -  °  (4,20) 

According  to  this  Lagrangian  criterion  of  stability,  the  least  stable 

root  is  that  root  which,  for  given  values  of  a   ,  U(y)   and  the  other 

R 

parameters,  has  the  highest  algebraic  value  of  y'  .   This  will  also  be 

R 

the  root  which  hes  the  highest  value  of  y   with  respect  to  the  original 

R 

fixed  axes . 

It  also  follows  from  the  foregoing  concept  that  different  streamlines 

exhibit  different  degrees  of  instability  (except  for  the  special  case 

where  a  =  0)  .   Of  particular  importance  is  the  least  stable  streamline, 
R 

that  is,  the  streamline  which,  for  given  values  of  y^,   and  aD  ,  has  the 

RL  R 

algebraically  largest  value  of  y'   .   Eq.  (4,20)  reveals  that  if  a   be 

RL  R 
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positive,  then  the  streamline  at  mid  channel  is  the  least  stable,  while 

if  a   be  negative,  then  the  streamline  at  the  wall  is  the  least  stable. 
R 

In  either  case  we  can  state  that  a  condition  of  incipient  instability  is 
reached  when 

▼k  =  °  (4-21) 

for  the  least  stable  streamline. 

From  the  foregoing  considerations  we  may  conclude  finally  that  the 
flow  is  stable,  or  at  least  neutrally  stable, 


when  aR  >  0  ,   if  Y^  <  -  f  aR 


when  a_  <   0  ,   if  yDT  <  0 


(A. 22) 


Notice  that  the  least  stable  root   Yr,T   ^n  Ecls«  (4. 22)  is  now  the 

KL 

ordinary  root  with  respect  to  stationary  axes;  the  moving  axes  have  served 

their  purpose  and  are  not  needed  further.   Also  notice  that  the  second  of 

Eqs.  (4.22)  is  identical  to  the  conventional  stability  criterion  of  Eq. 

(3.3).   Thus  only  if  a   >  0  do  we  need  to  modify  the  conventional 

R 

stability  criterion  in  the  manner  shown. 

The  generalized  Lagrangian  criterion  of  stability  as  developed  in 
this  section  and  summarized  concisely  in  Eqs.  (4.22),  forms  an  essential 
element  of  the  proposed  research.   This  generalization  of  the  stability 
criterion  represents  an  original  contribution  of  the  present  investigator 
and,  so  far  as  he  is  aware,  does  not  duplicate  the  published  work  of  any 
other  researcher. 

Moreover,  preliminary  results  calculated  by  this  method  and  reported 
by  Harrison,  Newby,  and  Johnston,  Ref s .  (6),  (7)  and  (8),  appear  reasonable 
and  are  encouraging.   They  show,  for  example,  that  the  introduction  of 
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negative  values  of  a   is  definitely  destabilizing  and  lowers  the  critical 

Reynolds  number  below  that  which  applies  for  the  conventional  case, 

a  =  0  .   This  is  a  significant  result.  Moreover,  it  is  a  result  that 

could,  not  have  been  discovered  if  the  investigation  had  been  guided  by 

the  overly  restrictive  criterion  of  neutral  stability  as  expressed  by 

Eqs.  (4.3). 

The  present  theory  also  shows  that  for  negative  values  of  a 

R 

instability  starts  at  the  walls  at  some  critical  Reynolds  number,  and  with 

gradually  increasing  Reynolds  number  it  spreads  progressively  inward. 

This  is  in  agreement  with  experimental  observation,  Schlicting,  Ref.  (3). 

Application  of  the  foregoing  concepts  permits  us  to  assign  the  four 

parameters  a  ,  a  ,  3_,  BT  arbitrarily  and  to  calculate  by  a  process  of 
R    I    R    I 

numerical  trial  and  correction  the  corresponding  Reynolds  number  for 
incipient  instability.  The  results  disclosed  by  calculations  of  this 
kind  for  a  suitable  range  of  values  of  the  above  four  parameters  can  be 
symbolized  in  the  form 

\  "  \  [V  V  br>  8i!  (4-23) 

Of  particular  interest  is  the  question  whether  there  exists  some 

particular  combination  of  parameters  a  ,  a  ,  3n,  3T  for  which  R 

K   I   R   I  e  . 

attains  an  absolute  minimum  value,  or  whether  a  ,  a  ,  3^  and  3T  > 

R    I    R         -L 

if  otherwise  unrestricted,  may  always  be  so  chosen  as  to  reduce  R 

below  any  pre-assigned  limit.   This  amounts  to  asking  whether,  or  in  what 
sense,  a  critical  Reynolds  number  may  be  said  to  exist  below  which  the  flow 
field  is  stable  to  all  three  dimensional  perturbations.   While  this  pro- 
blem has  been  solved  for  the  special  case  where  a  =  3^  ■  0  ,  it  remains 

R     R 

relatively  unexplored  for  the  far  more  complex  generalized  case. 
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5.   Transformation  of  Parameters  for  Plane  Poiseuille  Flow 

The  foregoing  development  has  shown  that  the  linearized  stability 
problem  of  plane  Poiseuille  flow  can  be  formulated  in  terms  of  the  para- 
meters a  ,  a  ,  6  ,  S  ,   and  R   .   In  this  section  we  introduce  an 

k   j.   k   i         e 

alternative  but  equivalent  set  of  parameters  which  are  here  symbolized  as 
A,  <J>,  9,  <  and  R  *  .   This  change  in  parameters  yields  an  important 

theoretical  result  and  also  permits  the  numerical  calculations  to  be 
significantly  reduced. 

The  fundamental  relations,  Eqs.  (2.11)  and  (2.12),  which  govern  the 
solution  for  the  eigenvalues  and  eigenf unctions  are  seen  to  contain  two 
important  complex  constants  as  shown  in  the  following  two  equations.   The 
left  side  of  each  equation  expresses  each  constant  in  terms  of  the  four 

original  parameters  a  ,  a  ,  8  ,  8T  and  the  right  side  in  terms  of  the 

K.    1    K    1 

four  new  parameters     A,    <f>,    9     and     <    .      Thus 


(CtR  +   ±   ai)2   +    (8R  +   i    V2    =   a2q   U 


(a_  +  i  aT)    =  KAe 

K  J. 


(5.1) 
(d>  +  9) 


Eqs.  (5.1)  suffice  to  fix  the  solution  for  the  four  new  parameters 

A,  <j>,  9  and  <  as  functions  of  the  four  original  parameters 

a  ,  a  ,  8^,  8T  or  vice  versa.   The  details  of  these  conversions  in 
R    X    R    1 

either  sense  are  given  by  Newby,   Ref.  (7),  and  are  not  repeated  here. 
Newby  also  gives  certain  auxiliary  relations  associated  with  this  trans- 
formation of  parameters  which  need  not  be  considered  in  the  present 
discussion. 
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It  is  useful,  however,  to  define  here  the  important  auxiliary 
parameter 

AR2  -  \     +  BR2  <5-2> 

The  relationships  symbolized  by  Eqs.  (5.1)  can  also  be  summarized 
graphically  in  the  complex  plane  as  shown  in  Fig.  5.1. 

The  vector  OP  in  the  diagram  represents  the  complex  quantity 

2  2 

a   .   The  vector  PQ  represents   8   •   The  geometry  of  the  diagram  is 

completely  fixed  when  the  coordinates  of  the  points  P  and  Q  are 
specified.   Notice  that  the  locations  of  these  two  points  suffice  to 
determine  the  four  characteristics  a  ,  a  ,  8-,  8T  or  A,  <J>,  Q     and 
k   which  fix  the  spatial  wave  form  of  the  perturbation. 

Consider  the  Reynolds  number  which  corresponds  to  incipient  insta- 
bility for  a  perturbation  whose  spatial  wave  form  has  been  specified. 
This  relation  is  symbolized  by  Eq.  (4.23)  which  we  repeat  here  for  easy 
reference,  namely, 

Rel  "  Rel  [V  V  6R>  SI]  (4-23) 

A  remarkable  and  important  feature  of  the  transformation  defined  by 
Eqs.  (5.1)  and  (5.2)  is  that  it  permits  the  relation  symbolized  by  Eq. 
(4.23)  to  be  simplified  significantly  as  explained  below. 

The  transformation  in  question  makes  use  of  a  special  reference 
case  which,  for  arbitrary  values  of  parameters  A,  <j>  and  0  ,   is 
characterized  by  a  particular  value  of  the  fourth  parameter,  namely, 

k  -»■  k*  =  1  (5.3) 
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It  can  also  be  shown  from  Eqs.  (5.1)  that  under  these  conditions  the 

2 
parameter  A    also  takes  on  a  corresponding  reference  value,  namely, 
R 


2     *2    2   i  2 

A   -*■  A    =  A   [|sin  Q|  +  cos  <j>J 


(5.4) 


For  the  above  reference  case,  the  solution  of  the  incipient  sta- 
bility problem  analogous  to  that  symbolized  by  Eq.  (4.23)  reduces  to  a 
form  which  may  be  symbolized  as  follows. 


R   +  R  *  =  R  *  [A,  <J>,  9] 
el    el     el 


(5.5) 


The  important  feature  of  Eq.  (5.5)  as  compared  with  Eq.  (4.23)  is 
that  the  number  of  independent  parameters  appearing  on  the  right  has  been 
reduced  from  four  to  three,  namely,  to  A,  $,   and  0  .   Owing  to  the 
length  and  complexity  of  the  matrix  calculations  involved,  this  reduc- 
tion from  four  to  three  independent  parameters  represents  a  very  great 
reduction  in  the  overall  computational  burden. 

The  problem  that  now  remains,  however,  is  how  to  generalize  from  the 
above  special  reference  case  to  the  more  general  case  in  which  parameters 

A,  <j>  and  0  remain  fixed  at  their  previous  values,  but  in  which  <   is 

2 
allowed  to  vary  arbitrarily.  Of  course  parameters  A„   and  R  _  also 

R        el 

vary  accordingly.   It  can  be  shown,  however,  that  these  variations  are 
governed  by  the  following  simple  analytical  relations,  namely, 


=  K 


-1 


(5.6) 
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The  detailed  derivation  of  the  relations  symbolized  by  Eqs.  (5.3) 
through  (5.6)  lies  outside  the  scope  of  this  discussion  but  is  fully 
documented  in  Ref.  (7). 

The  general  nature  of  the  relation  defined  by  Eq.  (5.6)  is  displayed 

graphically  in  Fig.  5.2.   Notice  that  for  fixed  values  of  A,  $  and  0  , 

2  -1 

an  increase  in  A    always  implies  a  corresponding  decrease  in  <  ' 
K 

2 
and  hence  a  decrease  in  R   .  Moreover,  if  A    be  allowed  to  increase 

e  J.  R 

Without  bound,  then  R    approaches  zero  in  the  limit.   On  the  other 

2 
hand  if  A    is  limited  to  some  specified  upper  bound,  then  R  _  cannot 
K  el 

drop  below  some  corresponding  lower  bound. 

The  foregoing  discussion  is  based  on  holding  parameters  A,  <£  and  0 

2 
fixed  and  considering  the  effects  of  possible  changes  in  A    .   Let  us 

R 

2 
next  consider  the  case  where  A    is  held  fixed  and  parameters  A,  <f> 

R 

and  0   are  free  to  vary.   It  follows  from  Eqs.  (5.5)  and  (5.6)  that 

there  now  exists  some  definite  combination  of  values  for  A,  <j>  and  0  , 

respectively,  for  which  R  _  attains  an  absolute  minimum  value.   We 

el 

term  this  the  critical  Reynolds  number  R  )    .   It  is  evident  that 

e  cr 

R  )    is  some  definite  function  of  the  single  remaining  independent 

2 
parameter,   A    .   One  important  and  specific  objective  of  the  proposed 

research  is  to  compute  this  function. 

2 
Eqs.  (5.5)  and  (5.6)  show  that  for  an  assigned  value  of  A    ,  the 

R 

critical  combination  of  values  of  A,  <j>  and  0  which  minimize  the  value 
of  R    cannot  be  found  analytically,  but  must  be  determined  by  a  syste- 
matic process  of  numerical  exploration  over  the  four  dimensional  domain 
whose  coordinates  are  A,  <j>,  0  and  R   .   Some  preliminary  numerical 
calculation  of  this  kind  has  already  been  accomplished  on  a  limited  scale, 
and  research  support  is  needed  for  the  purpose  of  extending  this  exploration 
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in  a  systematic  way.   This  task  forms  an  essential  part  of  determining 

2 
the  critical  Reynolds  number  as  a  definite  numerical  function  of  X 

R 

The  theory  developed  in  this  section  shows  that,  even  before  any 
specific  numerical  calculations  are  undertaken,  the  so-called  critical 

Reynolds  number  is  not  some  absolute  constant  as  it  is  commonly  conceived 

2 
to  be,  but  rather  some  definite  function  of  X        .   This  is  a  significant 

R 

new  concept  and  represents  an  original  contribution  by  the  present 
investigator. 

For  the  special  case  where 

AR2  =  0  (5.7) 

the  present  theory  and  transformations  reduce  to  the  conventional  and  well 
known  relations  associated  with  Squire's  theorem.   This  investigator  knows 
of  no  published  result,  however,  which  is  fully  equivalent  to  the  generalized 
form  of  these  transformation  as  outlined  in  this  section.   These  results 
are  believed  to  represent  an  original  contribution  by  the  present  investigator. 

The  analysis  of  the  generalized  criterion  of  stability  as  presented  in 
the  previous  section  has  shown  that  the  conventional  constraint  expressed  by 
Eq.  (5.7),  which  is  widely  regarded  as  a  necessary  condition  for  neutral  sta- 
bility, is  not  in  fact  necessary  at  all.   The  analysis  of  this  section  now 

2 
shows  that  positive  values  of  X  have  the  further  effect  of  decreasing 

R 

the  critical  Reynolds  number.   This  result  helps  explain  why  the  conventional 
analysis  tends  to  predict  values  of  critical  Reynolds  number  which  are  some- 
what too  high  and  how  slight  differences  in  experimental  set-up  or  test 

2 
conditions,  which  might  be  associated  with  different  values  of  X        ,  could 

K 

produce  apparent  inconsistencies  in  the  observed  critical  Reynolds  number0 
A  significant  point  is  that  this  explanation  of  discrepancies  is  still 
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based  on  a  relatively  straightforward  linearized  analysis  whereas  some 
investigators,  failing  to  find  any  simple  explanation  of  this  kind,  have 
felt  compelled  to  resort  to  non-linear  analyses,  with  all  of  the  enormous 
complexities  they  entail,  in  an  effort  to  make  sense  out  of  the  experi- 
mental observations. 
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6.   Stability  of  Pipe  Flow 

While  the  stability  of  plane  Poiseuille  flow  as  discussed  in  the 
previous  sections  of  this  report  is  of  fundamental  interest  because  of  the 
basic  simplicity  of  the  three  dimensional  rectilinear  geometry  involved, 
this  configuration  is  of  somewhat  limited  practical  importance,  and  it 
does  not  lend  itself  readily  to  experimental  study.   In  fact  the  amount  of 
experimental  data  available  for  this  case  is  extremely  limited. 

Pipe  flow  on  the  other  hand,  while  involving  a  geometry  which  is  only 
slightly  more  complex,  is  of  the  greatest  practical  importance.   Also, 
there  exists  an  extensive  body  of  experimental  data  on  the  stability  of 
pipe  flow.  Moreover,  the  discrepancy  between  the  conventional  theory  and 
experimental  observation  is  extreme  in  this  case  in  that  the  theory  fails 
completely  to  predict  the  experimentally  observed  instability.   These  are 
strong  reasons  for  generalizing  the  theoretical  analysis  for  pipe  flow  to 
include  three  dimensional  effects  in  a  manner  analogous  to  that  previously 
discussed  in  connection  with  plane  Poiseuille  flow.   Preliminary  calcula- 
tions based  on  such  a  generalized  analysis,  while  of  very  limited  scope, 
have  succeeded  for  the  first  time  in  disclosing  theoretical  instabilities 
in  pipe  flow.   This  is  a  noteworthy  accomplishment  and  warrants  a  sustained 
follow  up. 

Broadly  speaking,  the  pipe  flow  problem  can  be  analyzed  by  much  the 
same  method  previously  explained  in  connection  with  plane  Poiseuille  flow. 
Only  certain  details  are  different  in  the  two  cases.  An  extensive  analysis 
for  pipe  flow  is  given  by  Harrison,  Ref.  (6)  and  Johnston,  Ref.  (8),  and 
is  not  repeated  here.   However,  we  do  point  out  in  this  section  the  main 
respects  in  which  the  analysis  of  the  stability  of  pipe  flow  differs  from 
that  of  plane  Poiseuille  flow. 
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All  equations  are  nondimensionalized  by  utilizing  as  dimensional 
reference  parameters  the  density  p   of  the  fluid,  the  radius  a  of  the 
pipe  and  the  mean  volumetric  velocity  U  of  the  fluid. 

Cylindrical  coordinates  x,  r,  0  are  used  with  x  being  in  the 
direction  of  the  flow,  r  being  in  the  radial  direction  and  0  being  the 
angular  coordinate.   Dimensionless  coordinate  r  varies  from  zero  at  the 
axis  to  +1  at  the  pipe  wall.   Symbols  e  ,  e  ,  e   denote  unit  vectors 
in  the  three  coordinate  directions. 

It  is  well  known  that  the  dimensionless  velocity  U  of  the  mean  flow 
along  the  x  axis  is  given  by  the  simple  expression 

U  =  2(1  -  r2)  (6.1) 

It  is  of  interest  to  compare  this  relation  for  pipe  flow  with  the  corres- 
ponding relation  for  plane  flow  as  given  by  Eq.  (2.1).   Notice  that  the 
maximum  dimensionless  velocity  now  equals  2  as  compared  with  3/2  for  plane 
flow. 

The  vector  potential  of  the  perturbation  is  chosen  in  the  form 

+  +       +         +  (ax  +  i  n0  +  yt) 

W  =  [e  (o)  +  e  G  (r)  +  enH  (r)  ]  e  (6.2) 

n     x       r  n       0  n 

where  the  complex  constants  a  and  y  may  be  expanded  further  as  follows, 


a  =  a   +  i  a 


y  =  yR  +i  yx 


(6.3) 


Notice  that  the  coefficient  of   0  in  the  exponent  in  now  taken  not 
as  a  general  complex  quantity  3  ,  but  rather  as  the  purely  imaginary 
quantities  i  n  where  n  is  an  integer.   Thus 
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6  =  BR  +  i8=0  +  in 


3R  =  0  (6.4) 


3  =  n    n  =  0,  1,  2,  3,  ... 


The  reasons  for  restricting  the  complex  coefficient  8  in  this  way  is 
that  the  velocity  components,  and  hence  the  vector  potential  itself,  must 
be  strictly  periodic  with  respect  to  angle  0  .   This  is  a  fortunate  circum- 
stance since  it  eliminates  parameter   g   and  reduces   6T  to  an  integer, 
thus  simplifying  the  problem.   The  spatial  form  of  the  perturbation  is 

therefore  completely  defined  by  the  three  parameters  a  ,  a   and  n  . 

R   I 

As  Refs.  (6)  and  (8)  show  in  full  detail,  the  vorticity  transport 
equations  that  result  from  the  use  of  the  above  vector  potential  are  found 
to  reduce  to  the  following  form 

(D*Gn)  P\\  f^n)  (DGn 

[V      D4H        +    [M3'1      03H        +    [V>      D2H       +   IV 1      DH 

vn7  vny  vny  ^       n 


/  (d2g  }  (DG   ) 

\    2        (DZHJ         Ll    J   \mn\        LNo 


+  [M0-]  1M-t([H2']  P£j  +  [N^]  p»J+  [Ho-]  Q)  (6.5) 


where  symbol  D  now  denotes  the  differential  operator  d/dr  . 

The  eight  matrices   [Mf]   and  [Nf]  which  appear  in  Eq.  (6.5)  each  have 
two  rows  and  two  columns,  or  32  elements  in  all,  of  which  10  happen  to  be 
zeros.  All  elements  are  defined  in  detail  in  Refs.  (6)  and  (8).   These 
definitions  are  too  lengthy  to  be  repeated  here0   It  suffices  to  remark 
that  all  of  these  matrix  elements  are  known  complex  functions  of  the  three 
spatial  wave  form  parameters   a  ,  a   and  n  ,  of  Reynolds  number  R   , 
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and  of  the  radial  coordinate  r  .  Moreover,  the  coefficients  of  the 

3         4 
quantities  D  G   and  D  G   in  Eqs.  (6.5)  are  zeros,  so  that  these  equation 

are  of  only  second  order  in  G   whereas  they  are  of  fourth  order  in  H 

n  n 

Thus  only  two  boundary  conditions  are  needed  which  involve  G   or  its 

n 

derivatives,  while  four  are  needed  which  involve  H   or  its  derivatives. 

n 

An  interesting  feature  of  Eqs.  (6.5)  is  that,  in  general,  these  two 

equations  are  coupled,  that  is,  both  functions  G   and  H   and  their 

n        n 

various  derivatives  appear  in  both  equations.  Hence  the  equations  must 
be  solved  simultaneously.   For  the  special  case  n  =  0  ,  however,  these 
equations  uncouple.  One  of  them  becomes  an  eigenvalue  equation  of  fourth 
order  in  the  single  unknown  function  H   .   The  other  becomes  an  independent 
eigenvalue  equation  of  second  order  in  the  single  unknown  function  G   .   Th 
the  solution  for  n  =  0  breaks  down  into  two  independent  families  of  eigen- 
values and  eigenf unctions.   It  can  be  shown  that  one  of  these  families 
represents  perturbations  in  the  meridional  plane  9  =  constant  ,  the  other 
represents  perturbations  in  the  cross-sectional  plane  x  =  constant  .  On  tri 

other  hand  for  n  =  1,  2,  3  ...   the  eigenf unctions  G   and  H   which 

n        n 

correspond  to  a  particular  eigenvalue  remain  coupled  and  therefore  represent 

perturbations  of  a  truly  three  dimensional  character. 

Solution  of  the  basic  vorticity  transport  relations,  Eqs.  (6.5)  must 

be  found  subject  to  the  appropriate  boundary  conditions  at  the  wall.,   r  =  1 

and  at  the  axis,   r  =  0  .   The  boundary  conditions  are  discussed  in  some 

detail  in  the  next  section  of  this  report. 

If  numerical  values  of  a_,  a_  ,  n  and  R   be  specified,  the  real 

K   l  e 

part  Ydt   °f  t^e  least  stable  root  can  be  found  by  solution  of  the  vorticit 
transport  equations.   The  stability  can  then  be  determined  by  a  Lagrangian 
criterion  analogous  to  Eq.  (4.22).   In  other  words  the  flow  is  stable,  or 
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at  least  neutrally  stable, 


when  aR  >  0  ,   if  y^  <  -  2aR 
when  c^.<0  ,   if  yR  <  0 


Applications  of  the  foregoing  concepts  permits  us  to  assign  the  three 
parameters  a  ,  a   and  n  arbitrarily  and  to  calculate  by  a  process  of 
numerical  trial  and  correction  the  corresponding  Reynolds  number  for 
incipient  instability.   The  results  disclosed  by  calculation  of  this  kind 
for  a  suitable  range  of  values  of  the  above  three  parameters  can  be 
symbolized  in  the  form 

Rel  =  Rel  [V  V  nl  (6'7) 

Of  particular  interest  is  the  question  whether,  if  we  hold  a   con- 

R 

stant  at  some  arbitrary  value,  there  exists  some  particular  combination  of 

parameters  a_  and  n  for  which  R  _  attains  an  absolute  minimum  value. 
I  el 

We  term  this  the  critical  Reynolds  number  R  )    .   It  is  evident  that 

e  cr 

R  )  is  some  definite  function  of  the  single  remaining  independent  para- 
meter, a  .  One  important  and  specific  objective  of  the  proposed  research 
is  to  compute  this  function. 

While  the  stability  of  pipe  flow  has  been  extensively  investigated  for 

the  special  case  a  =  0  ,  it  nevertheless  remains  relatively  unexplored 

R 

for  the  more  general  case  in  which  a  ^  0  .  Moreover,  in  pipe  flow  certain 

R 

boundary  conditions  at  the  axis  become  so  problematical  as  to  raise  serious 
questions  about  the  validity  of  any  purported  solution  that  does  not  properly 
take  these  subtleties  into  account.   These  aspects  are  considered  in  the 
next  two  sections. 
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7.   Boundary  Conditions  for  Pipe  Flow- 

In  pipe  flow,  the  final  vorticity  transport  equations  are  of  second 

order  in  function  G   and  of  fourth  order  in  function  H   .   Hence  two 

n  n 

boundary  conditions  are  needed  with  respect  to  G   and  four  with  respect 

n 

to  H     It  turns  out  that  three  of  these  conditions  are  defined  at  the 
n 

wall,  r  =  1  ,  and  the  remaining  three  at  the  axis,  r  =  0  . 

The  boundary  conditions  at  the  wall,  r  =  1  ,  are  based  on  the  usual 
no  slip  restriction.   Thus,  for  all  values  of  index  n 

un(l)  «  0 


v  (1)  -  0 


w  (1)  =  0  J 

n 


(7.1) 


Inasmuch  as  the  governing  differential  equations  are  finally  expressed 

in  terms  of  the  two  basic  functions  G  (r)   and  H  (r)  ,  it  is  useful  to 

'  n  n 

rewrite  Eqs.  (7.1)  in  like  terms.   In  this  connection  recall  that  the 

-*-  •*- 

perturbation  velocity  v   is  related  to  the  vector  potential  W   as 
r  n  n 

follows. 


v  =  VXW  = 
n      n 


e  u  (r)  +  e  v  (r)  +  eA  w  (r) 

o     n 


x  n 


r  n 


(ax  +  inQ  +  yt) 
e  (7.2) 


From  Eqs.  (6.2)  and  (7.2)  we  readily  find  that  the  three  velocity  com- 
ponents are  given  by  the  following  expressions. 

u  =  DH  +  -  (H  -  in  G  ) 
n     n   r   n       n 


v  =  -  aH 
n       n 


(7.3) 


w  =  +  aG 
n        n 
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It  is  also  useful  in  the  present  context  to  expand  functions  G 

n 

and  H   in  Kaclaurin  series.   Thus 
n 

2  r2 

G     =  G   (o)    +  DG    (o)    r+DG(o)TT+   ... 
n  n  n  n  Z,  .      .. 

(7.4) 

2  2 

H     =  H   (o)    +  DH   (o)    r  +  D  H    (o)   ~r  +   . .  . 
n  n  n  n  2. , 


With  the  aid  of  Eqs.  (7.3)  and  (7.4)  it  becomes  a  straightforward 
task  to  rewrite  Eqs.  (7.1)  in  the  required  form.  The  following  results 
are  obtained  which  apply  for  all  values  of  index  n  ,  namely, 

Gn(l)  -  0 

Hn(l)  =  0  (7.5) 

DH  (1)  =  0 
n 

Derivation  of  the  remaining  three  boundary  conditions  which  apply  at 
the  axis  r  =  0  is  considerably  more  complicated.  We  start  by  observing 
that  Eq .  (7.2)  correctly  gives  the  perturbation  velocity  v   at  any  point 
which  is  not  on  the  axis  but  does  not  necessarily  define  a  meaningful 
velocity  at  a  point  on  the  axis.   To  express  the  velocity  on  the  axis  it  is 

not  sufficient  merely  to  set  r  =  0  in  Eq.  (7.3).   The  reason  is  that 

"*■   •*"         inO 
Eq.  (7.2)  contains  the  quantities   e  ,  eQ  and  e     which  are  functions 

of  9  ,  whereas  coordinate  0  ,  while  determinate  for  any  point  off  the 

axis,  is  indeterminate  at  the  axis  itself. 

We  may  circumvent  this  difficulty  by  starting  at  some  arbitrary  point 

near  but  not  on  the  axis,  and  proceeding  toward  the  axis  along  a  radial  line 

for  which  0  remains  constant.   In  the  limit  as   r  "*"  0  ,  the  final  location 

of  the  point  on  the  axis  is  independent  of  the  coordinate  0  of  the  line 
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along  which  that  point  is  approached.  Likewise,  the  velocity  v  at  that 
point  is  independent  of  the  coordinate  0  which  characterizes  the  path  of 
approach.   Hence  we  may  write 


30 


{r^O  ^n  (x'  r>  9'  t}}  =  °  (7>6) 


-)■ 


x 


In  applying  Eq.  (7.6)  to  Eq.  (7.2)  we  make  use  of  the  fact  that  e 

->       -> 
is  a  constant  but  e   and  en  are  functions  of  0  ,  so  that 

r      y 

3e 

k*r'  -  ° 

3e 

'  -  +  efi  (7.7) 


^30  / 

'^0 


,30  /  =  "  Gr 


It  now  follows  easily  from  Eqs .  (7.2),  (7.6)  and  (7.7)  that  for 
arbitrary  0  and  arbitrary  n 

s 

n  u  (o)  =  0 
n 

i  n  v  (o)  -  w  (o)  =  0  (7.8) 

n       n 

v  (o)  +  in  w  (o)  =  0 
n         n 

Eqs.  (7.8)  stipulate  the  conditions  that  must  be  satisfied  in  order 
that  Eq.  (7.2)  shall  define  a  unique  and  determinate  value  of  v   in  the 
limit  as  r  "*"  0  ,   The  value  of  0  at  this  limit  may  be  chosen  arbitrarily 
without  affecting  the  resulting  value  of  v   . 

Eqs.  (7.8)  break  down  into  three  distinct  cases,  namely,   n  =  0  , 
n  =  1   and   n  ;>  2  .   Thus 
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For  n  =  0 


v  (o)  =  0 
o 


w  (o)  =  0 
o 


(7.9) 


Notice  that  u  (o)   does  not  vanish  in  this  case 
o 


For  n  =  1 


u1(o)  =  0 


V;L(o)  +  i  w1(o)  =  0 


(7.10) 


Notice  that  v  (o)   and  w  (o)  ,  although  mutually  constrained,  do 

not  vanish  in  this  case. 

For  n  =  2,  3,  4,  5,  ... 

u  (o)  =  0 
n 

v  (o)  =  0  ,  (7.11) 

n  ' 

w  (o)  =  0 
n 


It  is  significant  that  u  (o)  ,  v  (o)   and  w  (o)   are  the  only  non- 
zero velocity  components  at  the  axis.   Thus  the  cases  n  =  0   and  n  =  1 
suffice  to  account  fully  for  the  three  non  zero  perturbation  velocity  com- 
ponents that  can  occur  at  this  location. 

It  is  now  useful  to  express  Eqs.  (7.9),  (7,10)  and  (7,11)  in  terms  of 

functions  G   and  H   .   The  following  results  are  obtained 

n        n  ° 


case. 


For  n  =  0 

G  (o)  =  0 
°  (7.12) 

H  (o)  =  0 
o 

Note  that  one  additional  boundary  condition  is  still  needed  for  this 
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For  n  =  1 


H^o)  -  i  G1(o)  =  0 


2DH  (o)  -  i  DG  (o)  =  0 


(7.13) 


Note  that  one  additional  boundary  condition  is  still  needed  for  this 

tease  also. 

Notice  in  addition  that  this  case  is  exceptional  in  that  it  represents 

(the  only  value  of  n  for  which  H  (o)   and  G  (o)   do  not  vanish, 

n  n       

For  n=2,  3,  4,  5,  ... 
Gn(o)  =  o 

H  (o)  =  0  (7.14) 

n 

2DH  (o)  -  i  n  DG  (o)  =  0 
n  n 

It  is  now  clear  that  for  n  ^_  2    ,  Eqs.  (7.5)  and  (7.14)  fix  all  six 
of  the  boundary  conditions  needed  to  define  a  determinate  solution.   On  the 
other  hand  the  cases  n  =  o  and  n  =  1  each  still  lack  one  final  necessary 
boundary  condition. 

The  one  remaining  boundary  condition  needed  for  the  case  n  =  0  can 

readily  be  deduced  from  well  known  physical  laws,   Consider  a  small  fluid 

element  of  cylindrical  shape,  of  length  Sx  and  radius   5r  ,  whose  axis 

coincides  with  the  x  coordinate  axis.   Let  x    denote  the  longitudinal 

xr  to 

shear  stress  acting  on  the  curved  surface  of  the  element,  let  p  be  the 

TVii 

pressure,   p   the  density  and   ("77")   the  axial  acceleration.  We  may  write 
the  equation  of  motion  in  the  x  direction  for  this  element  as  follows, 

T   27r5r  5x  -  (fj)  5xtt  (6r)2  =  p  Tr(5r)2  6x  (~)  (7.15) 

xl  ox  at 
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Proceeding  to  the  limit   5r  "*"  0  ,  we  observe  that  the  last  two  terms 
of  Eq.  (7.15)  are  of  higher  order  than  the  first.  Hence 

lim  x   =  lim  x   =  0  ,_..-. 

xr       xr  (7.16) 

5r  *  0    r  ■*■  0 

It  then  follows  from  the  known  relation  between  stress  and  strain 
rate  for  a  viscous  fluid  that  the  corresponding  shear  strain  rate  also 
vanishes  along  the  axis.   It  can  readily  be  shown  that  this  requires  in 
turn  that 


lim^Du  +  ccv}  .  0  (7_17) 


Upon  expressing  this  in  terms  of  functions  G   and  H   and  upon 
making  use  of  Eqs.  (7.9),  we  find  that  the  vanishing  of  the  shear  strain 
rate  requires  finally  that 

D2H  (o)  =  0  (7.18) 

o 

This  fixes  the  last  boundary  condition  needed  for  the  case  n  =  0  . 

All  of  the  boundary  conditions  required  for  the  pipe  flow  problem 
are  therefore  determined  above  except  one.   This  is  the  last  boundary  con- 
dition needed  to  define  a  determinate  solution  for  the  important  case 
n  =  1  . 

Unfortunately,  this  last  required  boundary  condition  is  also  the  most 

difficult  to  formulate  in  a  clear,  unequivocal  and  rigorous  manner.   The 

axis  represents  a  singular  line  along  which  many  of  the  usual  analytical 

expressions  for  quantities  like  v,  w,  Vxw  and  so  on  break  down.   The 

functions  G   and  H   and  their  various  derivatives,  especially  the 
n        n 

higher  derivatives,  are  not  necessarily  continuous  at  r  =  0  ,   Conse- 
quently, intuitive  notions  based  on  various  symmetry  relations  and 
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continuity  assumptions  become  unreliable  guides. 

We  call  attention  at  this  point  to  the  analogy  that  exists  between 
perturbation  vorticity  oj   and  perturbation  velocity  v  in  this  flow  field. 
Thus 

(ax  +  inO  +  yt) 


"n  =  V  X  \   =  lex  5(r)  +  er  n(r)  +  ^0  c(r)^  e  (7.19) 

The  analogy  between  this  relation  and  Eq.  (7.2)  implies  that  if 
Eq.  (7.19)  is  to  remain  valid  at  r  =  0  ,  then  we  must  have,  by  analogy 
with  Eqs.  (7.8) , 

n  Cn(o)  =  0 

i  n  nn(o)  -  cn(o)  =  0  (7>2Q) 

n  (o)  +  in  c  (o)  =  0 

n         n 

In  particular  for  n  =  1  ,  these  reduce  to 


51(o)  =  0 


n1(o)  +  i  C1(o)  =  0 


(7.21) 


Of  course  Eqs.  (7.13)  already  supply  two  of  the  three  needed  boundary 
conditions,  and  Eqs.  (7.21)  are  under  consideration  here  only  in  connection 
with  the  one  additional  constraint  that  is  still  lacking.   If  we  invoke  the 
first  of  Eqs.  (7.21)  for  this  purpose,  the  following  constraint  is  obtained, 
namely, 


DH  (o)  -  i  2DG  (o)  =0  (7.22) 
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On  the  other  hand,  if  we  invoke  the  second  of  Eqs.  (7.21)  for  the  same 
purpose,  a  quite  different  result  is  obtained,  namely 

D^Co)  -  i  D2G1(o)  =  0  (7.23) 

Unfortunately,  with  only  one  degree  of  freedom  still  at  our  disposal, 
there  is  no  way  to  satisfy  both  Eq.  (7.22)  and  Eq.  (7.23)  simultaneously. 
Hence  both  of  Eqs.  (7.21)  cannot  be  satisfied  simultaneously  either.   This 
simply  means  that  for  the  case  in  question,   n  =  1  ,  Eq.  (7.19)  inevitably 
breaks  down  at  r  =  0  .   Moreover,  the  attempt  to  utilize  the  analogy  between 
Eq.  (7.19)  and  (7.2)  leads  to  an  overspecification  of  the  needed  final 
boundary  condition. 

Further  study  is  needed  in  order  to  reformulate  the  required  final 
boundary  condition  at  the  axis  in  a  way  that  is  unambiguous  and  free  of 
arbitrary  assumptions.   This  is  an  important  question,  and  the  difficulties 
involved  in  correctly  formulating  this  particular  boundary  condition  might 
perhaps  account  for  the  failure  of  the  conventional  theory  to  disclose  the 
known  instability  of  pipe  flow.   A  cricial  review  of  this  question  will  form 
an  important  part  of  the  proposed  research. 

For  the  present,  we  shall  tentatively  assume  that  the  first  of  Eqs.  (7.21 
provides  the  appropriate  final  boundary  condition,  and  that  the  second  of 
these  equations  may  be  discarded.   This  amounts  to  retaining  Eq.  (7.22)  and 
discarding  Eq.  (7.23). 

Upon  combining  Eqs.  (7.13)  and  (7.22),  we  find  that  the  three  final 
boundary  conditions  at  the  axis  simplify  for  this  case  to  the  following  form 
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For  n  =  1 


DG  (o)  =  0 


H1(o)  -  i  G1(o)  =  0  (7.24) 


DH  (o)  =  0 


The  final  overall  boundary  conditions  for  all  values  of  n  are  now 
as  summarized  in  Eqs.  (7.25)  through  (7.28)  of  Table  7.1. 
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Table  7,1  Boundary  Conditions  for  Pipe  Flow 


For  n  =  1 


G  (1)  =  0 

o 


G  (o)  =  0 
o 


(7.25) 


H  (1)  -  0 


DH  (1)  =  0 
o 


H  (o)  =  0 
o 


D  H  (o)  =  e 
o 


(7.26) 


This  case  is  exceptional  in  that  the  solutions  for  functions  G 

i 

and  H   become  uncoupled  for  n  =  0  , 
n 

For  n  =  1 


Gx(l)  =  0 


DG  (o)  =  0 


H^l)  =  0 


H^o)  -  i  G1(o)  =  0 


(7.27) 


DH  (1)  =  0 


DH1(o)  =  0 


This  case  is  exceptional  in  that  the  quantities  G  (o)   and  H  (o) 

n  n 

become  non-zero   for     n  =  1    . 


For     n  =   2,    3,    4,    5, 


Gn(l)    =   0 

G  (o)  =  0 
n 

Hn(l)    -   0 

H  (o)  =  0 
n 

DH    (1)    =   0 
n 

2DH  (o)  -  i 
n 

(7.28) 


n 
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3.   Radial  Symmetry  Relations  and  Continuity  Conditions  at  Axis 

It  is  quite  common  in  dealing  with  a  variety  of  physical  problems 
tfhich  involve  a  cylindrical  geometry  to  derive  certain  boundary  conditions 
at  the  axis  from  considerations  of  symmetry.  Hence  it  is  natural  to  in- 
quire whether  any  such  constraints  apply  also  to  the  pipe  flow  problem. 
The  main  purpose  of  the  discussion  in  this  section  is  to  show  that  although 
the  solution  of  the  present  problem  does  indeed  satisfy  certain  relations 
of  radial  symmetry,  the  attempt  to  derive  the  boundary  conditions  on  this 
basis  involves  certain  hidden  pitfalls.  We  may  even  conjecture  that  such 
pitfalls  might  be  a  factor  in  the  failure  of  conventional  efforts  to  find 
any  theoretical  instability  in  pipe  flow,  although  we  know  from  experience 
that  such  instability  does  in  fact  occur. 

With  reference  to  symmetry,  consider  an  arbitrary  point  P  with 
coordinates  x,  r,  0  .   A  diametrically  opposite  point  Pf   may  now  be 
located  with  coordinates  x,  r,  (0  +  ir)  .  While  negative  values  of  coordi- 
nate r  are  not  ordinarily  used,  we  can  elect  in  this  instance  to  employ 
them,  and  to  designate  the  coordinates  of  point  P'   in  the  alternative 
but  equivalent  form  x,  -r,  0  .   Of  course  the  perturbation  velocity  at 
point  P'  ,  call  it  v   ,  must  be  the  same  regardless  of  which  of  the  above 
two  alternative  coordinate  descriptions  are  employed.   Thus 

v   [x,  -r,  0,  t]  =  v   [x,  +r,  (0  +  TT),t]  (8.1) 

Harrison,  Ref.  (6),  has  shown  that  as  a  consequence  of  the  relation 
symbolized  by  Eq.  (8.1),  the  functions   G   and  H   must  satisfy  symmetry 
constraints  of  the  form 
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G  („r)  =  -  (~l)n  Gt+r) 
n  n 

(8.2) 

H  (-r)  =  -  (-l)n  H  (+r) 
n  n 

This  tells  us  that  if  n  be  even,  functions  G   and  H   have  odd 

n        n 

symmetry  with  respect  to  r  ,  while  if  n  be  odd,  G   and  H   have  even 

n        n 

symmetry. 

Consequently,  for  n  even, 

G  (-r)  =  -  G(+r) 
n         n 

D2G  (-r)  =  -  D2G  (+r)  (8.3) 

n  n 

D4G  (-r)  =  -  D4G  (+r)    etc. 
n  n 

and  similarly  for  function  H   and  its  even  derivatives. 

n 

Likewise,  for  n  odd, 

DG  (-r)  =  -  DG  (+r) 

(8.4) 

D3G  (-r)  =  -  D3G  (+r)    etc. 
n  n 

and  similarly  for  the  odd  derivatives  of  function  H 

n 

If  we  now  assume  that  the  functions  and  derivatives  in  question  are 
also  continuous  at  r  =  0  ,  then  the  above  relations  have  the  following 
further  consequences.- 

For  n   even 

G  (o)  =  0  H  (o)  =  0 

n  n 

D2G  (o)  =  0  D2H  (o)  =  0  (8.5) 

n  n 

D4G  (o)  =  0  D4H  (o)  =  0 

n  n 

etc.  etc. 
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For  n  odd 

DG  (o)  =0  DH  (o)  =  0 

(8.6) 

D3G  (o)  =  0  D3H  (o)  =  0 

n  n 

etc.  etc. 

An  error  that  is  sometimes  made  is  to  assume  that  the  known  symmetry 
relations  expressed  by  Eqs.  (8.2)  necessarily  imply  the  consequences 
summarized  in  Eqs.  (8.5)  and  (8.6).   In  fact,  however,  these  consequences 
apply  if  and  only  if  the  quantities  in  question  are  known  to  be  continuous 
at  r  =  0  .   Inasmuch  as  the  mathematical  functions  ordinarily  encountered 
in  connection  with  incompressible  flows  are  usually  highly  continuous,  it 
is  easy  to  overlook  the  fact  that  Eqs.  (8.5)  and  (8.6)  actually  involve  a 
hidden  assumption  of  continuity  that  may  or  may  not  be  justified  in  the 
present  instance. 

The  situation  is  further  clouded  by  the  fact  that  a  number  of  the 
boundary  conditions  derived  in  the  previous  section  from  entirely  different 
considerations,  and  summarized  in  Table  7.1,  happen  to  coincide  with  the 
symmetry/ continuity  relations  of  Eqs.  (8.5)  and  (8.6).   These  are  the 
following 

For  n  =  0 

G  (o)  =  0 
o 

H  (o)  =  0  (8.7) 

o 

D2H  (o)  =  0 
o 
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For  n  =  1 


DG  (o)  =  0 


DH  (o)  =  0 


(8.8) 


For  n  =  2,  4,  6,  ... 


G  (o)  =  0 
n 


H  (o)  =  0 
n 


(8.9) 


Moreover,  for  the  case  n=3,  5,  7,  ...  ,  if  Eqs.  (8.6)  could  be 
assumed  to  apply,  they  would  also  ensure  that  the  condition 

2DH  (o)  -  in  DG  (o)  =  0  (8.10) 

n  n 

is  satisfied  identically. 

The  correct  interpretation  of  Eqs.  (8.7),  (8.8)  and  (8.9)  would  seem 
to  be  that  these  relations,  derived  earlier  on  quite  independent  grounds, 
ensure  that  the  corresponding  particular  quantities  must  also  be  continuous 
at  r  =  0  .   But  this  restricted  conclusion  does  not  warrant  the  indiscri- 
minate extension  of  Eqs.  (8.5)  and  (8.6)  to  other  cases  without  further 
justification. 

Furthermore,  there  are  a  number  of  boundary  conditions  listed  in 
Table  7.1  which,  while  not  necessarily  incompatible  with  Eqs,  (8.5)  and 
(8.6),  are  entirely  unrelated  to  the  latter.   These  are  the  following 

For  n  =  1 

H  (o)  -  i  G1(o)  =  0  (8.10) 

For  n  =  2,  4,  6,  ... 

2DH  (o)  -  i  n  DG  (o)  =  0  (8.11) 

n  n 
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For  n  =  1 ,  3 ,  5 ,  , ,  . 

G  (o)  =  0 
11  (8.12) 

H  (o)  =  0 

n 

If  we  proposed  for  some  reason  to  derive  the  boundary  conditions 
exclusively  from  symmetry  relations  and  continuity  assumptions,  it  would 
be  necessary  to  abandon  the  conditions  listed  in  Eqs.  (8.10),  (8.11)  and 
(8.12),  and  to  substitute  for  them  other  relations  drawn  from  Eqs.  (8.5) 
(8.6).   At  the  present  time  no  adequate  physical  or  mathematical  justifi- 
cation can  be  offered  for  such  a  change. 

In  an  effort  to  clarify  some  of  the  conceptual  complications 
associated  with  Eqs,  (8.5)  and  (8,6),  Johnston  in  Ref.  (8)  investigated 
the  constraints  that  would  theoretically  have  to  be  satisfied  in  order  to 
permit  the  direct  substitution  r  =  0  into  the  vorticity  transport  equa- 
tions.  His  principal  result  pertains  to  the  case  n  _>,  2  .   He  found  for 

this  case  that  functions  G   and  H   and  their  first  four  derivatives 

n        n 

would  all  have  to  vanish  at  r  =  0  to  allow  this.   His  result  involves 
ten  constraints  in  all.   Unfortunately,  with  only  three  degrees  of  freedom 
available  at  the  axis,  only  three  of  these  ten  constraints  can  actually  be 
enforced.   While  the  details  are  different  for  n  =  o  and  n  =  1  ,  the 
number  of  constraints  that  would  theoretically  be  required  for  this  pur- 
pose always  exceeds  the  three  available  degrees  of  freedom.   Hence  we 
conclude  that  the  simple  substitution  r  =  0  into  the  vorticity  transport 
equations,  in  their  original  analytical  form,  is  never  permissible.   In 
other  words,  the  axis  remains  irreducibly  a  singular  point,  and  Johnston's 
analysis,  while  correct  in  itself,  does  not  provide  the  three  and  only 
three  boundary  conditions  needed  at  that  point. 
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In  view  of  the  various  considerations  discussed  above,  it  is 
believed  that  Table  7.1  provides  the  best  overall  formulation  of  the 
boundary  conditions  that  is  available  at  this  time. 
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9.   Status  of  Pipe  Flow  Stability  Calculations 

Harrison,  Ref.  (6),  correctly  formulated  the  basic  equations  for 
pipe  flow  and  the  boundary  conditions  for  n  =  0  and  n  =  1  .  Unfortu- 
nately, he  was  unable  to  compute  results  for  these  cases  owing  to  a  minor 
error  in  his  computer  program.   Johnston,  Ref.  (8),  corrected  this  error 
and  computed  valid  results  for  n  =  0  ,  but  the  scope  of  his  calculations 
was  very  limited.   Johnston,  in  turn,  formulated  the  boundary  conditions 
incorrectly  for  all  cases  other  than  n  =  0  .   Hence  a  specific  aim  of 
the  proposed  research  is  to  amend  and  extend  the  above  calculations, 
using  the  corrected  boundary  conditions  of  Table  7.1.   It  is  felt  that 
the  revised  boundary  conditions  for  n  =  1  are  particularly  likely  to 
disclose  instabilities  that  have  not  appeared  before. 

Despite  the  error  in  boundary  conditions  involved  for  some  of  the 

cases,  Johnston's  results  in  Ref.  (8)  suffice  to  indicate  the  strong 

probability  that  small  negative  values  of  a   are  definitely  destabilizing, 

R 

An  important  aim  of  the  proposed  research  is  to  explore  this  question  more 
fully  and  adequately.   An  even  more  exciting  possibility  is  that  with  the 
revised  boundary  conditions  discussed  above  for   n  =  1  ,  instability  might 

possibly  be  found  even  for  a  =  0  .   If  so,  this  result  would  represent 

R 

a  very  significant  achievement. 
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10.   Summary  of  Research  Objectives 

The  first  general  objective  of  the  research  outlined  in  this  report 
is  to  carry  out  the  calculations  described  at  least  to  the  extent  necessary 
to  provide  significant  support  for  the  overall  validity  of  the  present 
theory  and  for  the  various  innovations  that  it  embodies.   For  the  case  of 
plane  flow  this  entails  exploring  some  typical  portion  of  the  boundary  of 
incipient  instability  as  expressed  symbolically  in  the  form 

Rei*  =  Rei*  [A'  *'  0]  t10'15 

For  the  case  of  pipe  flow,  it  entails  a  corresponding  exploration  of  the 
form 

Rel  =  Rel  [V  V  n]  (10*2) 

A  particular  objective  of  the  work  on  pipe  flow  is  to  confirm  the 
adequacy  of  the  revised  boundary  conditions  at  the  axis  as  summarized  in 
Table  7.1.   The  case  n  =  1  warrants  special  attention  as  having  possibly 
the  greatest  potential  for  disclosing  instabilities  not  detected  previously, 

The  ultimate  objective  is  to  carry  the  above  calculations  forward  to 
the  extent  that  it  becomes  possible  to  plot  the  two  important  functions 
symbolized  below  for  plane  flow  and  pipe  flow,  respectively,  namely, 

R     =  R     a  )    plane  flow  (10.3) 

e       e      R 
cr      cr 

R     =  R     (a  )    pipe  flow  (10.4) 

G  Q.  L\ 

cr      cr 

The  demonstration  of  the  very  existance  of  these  two  functions  would 
represent  a  noteworthy  and  fundamental  new  contribution  to  the  theory  of 
hydrodynamic  stability. 
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An  essential  collateral  objective  at  all  stages  of  the  research  is 
to  report  progress  periodically  in  the  technical  literature.   Enough  pro- 
gress has  already  been  achieved  to  warrant  at  least  one  initial  technical 
paper  of  this  kind.   Support  is  needed  as  much  for  preparing  such  a  paper 
on  progress  already  made  as  for  continuing  further  analysis  and  calculations, 
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